Viscosity of High Energy Nuclear Fluids 



V. Parihar 

Physics Department, Boston University, 



Boston MA 02215 



A.Widom, D. Drosdoff 

Physics Department, Northeastern University, Boston MA 02115 



o 
o 

(N 



in 



> 

ON 

o 
m 
O 

o 



o 

=5 



X 



Y.N. Srivastava 
Physics Department & INFN, University of Perugia, Perugia IT 

Relativistic high energy heavy ion collision cross sections have been interpreted in terms of almost 
ideal liquid droplets of nuclear matter. The experimental low viscosity of these nuclear fluids have 
been of considerable recent quantum chromodynamic interest. The viscosity is here discussed in 
terms of the string fragmentation models wherein the temperature dependence of the nuclear fluid 
viscosity obeys the Vogel-Fulcher-Tammann law. 
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I. INTRODUCTION 

There has been considerable interest in probing ex- 
perimentally a theoretical quark-gluon plasma phase 
transition by employing ultra-relativistic two heavy nu- 
clei scattering. It is not presently clear whether or 
not a deconfined plasma has actually been observed. 
However, such scattering experiments have apparently 
produced nuclear matter at very high energy density 
wherein the matter acts as a relativistic almost ideal 
fluidp], 0, H, 0b i.e. a quark-gluon matter fluid with very 
low viscosity]! [E 0, H ■ 

In previous work[^, we argued from the string 
fragmentation model, that a high energy density fluid 
may exhibit a glass-like behavior with meson and baryon 
strings playing a role analogous to polymers in some con- 
densed matter glasses. At very high temperatures, the 
viscosity of the stringy fluid is actually quite small. It 
is this regime (of an almost ideal fluid) that is of central 
importance in ultra-relativistic heavy ion scattering. Our 
purpose is to study how the viscosity of resulting nuclear 
matter varies with energy density, pressure and baryon 
density. 

In Sec|TT]we review the Kubo formula for the frequency 
dependent viscosity of a fluid and derive sum rules which 
allow for the determination of the zero frequency trans- 
verse viscosity in terms of the relaxation time for the 
non-diagonal elements of the pressure tensor. In Sec lIIIl 
we consider the quantized unit of circulation which con- 
tributes to turbulent eddy currents in the fluid. A bound 
for the Reynolds number of quantum turbulence leads 
to a quantum lower bound for the viscosity of the type 
that has been previously coniecturedpd|. We review the 
low energy well established liquid drop model in Sec lIVI 
and argue that the low energy density nuclear fluid is 
highly viscous. We turn to high energy nuclear matter 
in Sec|V] and in particular discuss the formal conformal 
symmetry associated with glue. This system is thought 
to describe high energy QCD nuclear matter in the limit 



of low baryon density. In Sec lVIl we review the notion of 
the QCD inspired string model and employ in Sec lVIII the 
string fragment ion picture to compute the viscosity of 
high temperature nuclear matter. The viscosity reaches 
small values in the ultra high temperature limit. In the 
concluding Sec lVIIIl we compare the string fragmenta- 
tion (jet) picture with the quark gluon plasma view of 
high temperature nuclear matter. 



II. PRESSURE FLUCTUATIONS 

For an isotropic fluid in thermal equilibrium, the off 
diagonal elements of the pressure tensor (say P xy (r,t)) 
fluctuations are described by the correlation function, 
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K{v-v',t) = -J (P xy (r,t)P xy (r',~i\))d\, (1) 

wherein (3 = (H/ksT). The complex frequency depen- 
dent fluid shear viscosity 77(C ) is determined by the Kubo 
formalism!!!, d Q M, UM 



V(0 



Q{t)= J K(r,t)d 3 r, 
g{t)e Kt dt, 3m C > 0. 



(2) 



With p as the mass density of the fluid and Coo as the 
velocity of a finite high frequency transverse sound, one 
finds the dispersion relations 
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together with the sum rule 
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In the low frequency limit 

ri = lim ner)(oj + i0+) = I Q(t)dt, (5) 



the relaxation time r for the decay of the correlation 
function 
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uniquely determines the shear viscosity via[l7L fl8l IT9l. l20l| 



V = P C oo T - 



(7) 



To find the viscosity at zero frequency for a fluid of mass 
density p employing the Kubo formalism, one must deter- 
mine the velocity of transverse sound Coo and the relax- 
ation time r of the off diagonal elements of the pressure. 



III. QUANTUM TURBULENCE 

The following arguments, first applied by Feynman [2l| 
to superfluid Helium, hold true for any fluid whose con- 
stituent particles all have mass M. For flows with eddy 
currents on a velocity scale V and a length scale L, the 
circulation of the eddy velocity field is given by 



v • dr = VL 



(8) 



From quantum fluid mechanics, the Feynman number T 
of quantum vortices in a circulating eddy is given by the 
Bohr quantization of circulation 



M (tvdr= ipdr = 2nhJ r ; 



i.e. the Feynman number is 



MVL 
2nh 



(9) 



(10) 



On the other hand, the Reynolds number for a circulating 
eddy is well known to be 
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pVL 



fill 



yielding a ratio which depends only on intrinsic fluid 
properties and not on the nature of the eddy flow; 



T Mr] 



K 2irfip 



(12) 



If p = Mn, with n as the number of constituent particles 
per unit volume, then 



T T) 
1Z 2nhn 



(non — relativistic). 



(13) 



Since the Feynman number of quantized vortices in an 
eddy is larger than the Reynolds number of the eddy, we 
have the quantum viscosity inequality 



T] > 2Tthn. 



(14) 



The above lower bound is found to apply to all known 
chemically pure substances. For example, for pure water 
T w 45ft > TZ. 



IV. OLD LIQUID DROP MODEL 

The first nuclear fluid model was employed in a non- 
relativistic context. In the liquid drop model, one has a 
momentum sphere of filled quasi-particle states out to a 
Fermi momentum pp = hkp. The number of nucleons 
A in a spherical nucleus of radius R can be computed 
by filling up a Fermi sphere in momentum space and a 
sphere in position space 



A = 



4 /47rfc| 



4vri? 3 



8(k F Rf 
9?r 



(15) 



The liquid drop nucleus is nearly incompressible so that 
the radius obeys 



R = aA 1 / 3 wherein 



1.2 x 10 _13 cm. 



(16) 



The Fermi velocity for a nearly ideal Fermi fluid is given 
by 



1 dEr, 



h dkf 



!)7ry /8 _R_ 
Ma 



(17) 



Note that (vp/c) w 0.27 which implies that the nucleons 
within the droplet are not entirely non-relativistic. 

The density of energy states per unit volume at the 
Fermi surface is given by 



9F 



(2tt) 3 



F dE F J 



) = A (3l) 

7T 2 \tiVF J 



(18) 



The specific heat per unit volume of a Fermi liquid is 
given by 



T- 



, ds 
dT 



ir 2 k 2 B g F 



T 



as T -> 0. (19) 



Which implies a low temperature entropy per unit vol- 
ume 



ir 2 k 2 B g F 



T = 



n k * 9F ) * (20) 
3 J ds 



wherein e is the excitation energy per unit volume. For 
weakly excited states the entropy per unit volume obeys 



2ir 2 g F 



1/2 



Ve = 2k B k F 



3Hv F 



(21) 
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In terms of the total droplet entropy and energy, the 
thermodynamic equation of state for excitation energy 
E weakly above the ground state energy E a is given by 



S(E) _ i A 
kn 3 



Trk 2 F a 3 ^ 1/2 
Uvf 



E-E 0A 



A 



(22) 



The density of excited energy levels for the droplet as a 
whole then reads 



D(E) w D(E 0A ) exp 



S(E) 



(23) 



which is experimentally reasonably accurate for heavy 
nuclei; a logarithmic plot of the density of excited states 
above the ground state yields a square root behavior [22l|. 
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D(E) 



D(E 0A ) 



IE-E 0A 



E A 



as E^ E aA + 0+. (24) 



In the opposite limit E > E 0A , one finds S(E) w E/T 
wherein To is the Hagedorn temperature [23|. 

The viscosity of the Fermi fluid in the liquid drop is 
quite high. This may be understood by writing Eq. ([7]) in 
the form 



V 



VpT c . 



(25) 



wherein r c is an appropriate collision time for quasi- 
particle scattering. From two body scattering phase 
space considerations, one finds that t c cx T -2 which im- 
plies a low nuclear energy liquid drop viscosity given by 



(26) 



in which K is a dimensionless quantity of order unity, 
depending on the angular variation of the two body cross 
section. At very low temperatures, the viscosity grows 
very high. It will be discussed in what follows that as the 
temperature grows ever larger, the viscosity diminishes to 
smaller values until it is quantum limited and the fluid is 
nearly ideal. 



V. GLUE AND CONFORMAL SYMMETRY 

If G a ^ v denotes the gluon field, 

G% = d,A a u - d v A« + (£) tfXAl, (27) 
then the action describing pure glue reads 

Under the conformal scale changes in space time 

x —>■ x = Ax, 

a;(x) - A;(x) = jA«(x), 

G«A X ) - = ±G%{x), (29) 



the action is left invariant; i.e. W — > W — W. Such 
symmetry dictates that the trace of the energy-pressure 
tensor 



-e = T^ = 'SP-u 



(30) 



is null. In Eq. (|30p . u represents the energy per unit vol- 
ume and P represents the pressure. Formally, the energy- 
pressure tensor for glue is given by 

Tpv"* — ^5 A<T G° (T G ai/ A — ^ff^i/G aM<T G QAIcr ^ , (31) 

which implies = and which formally verifies the con- 
formal invariance properties 



e 



Glue 



Glue 



= o. 



(32) 



On the other hand, for a theory with massive quarks the 
trace of the energy-pressure tensor is 



J2m f c 2 q f q f . 
f 



(33) 



In the presence of Gluons, the trace of the energy- 
pressure tensor can reassert its strength via 
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E 



m k c q k q k 



(34) 



For example, for a heavy quark, say Q with mass M, in 
the presence of a low energy gluon field configuration one 
finds the trace anomaly [24| 



9 Q - Mc 2 (0; G \QQ\ 0; G) = - (j^J G^G*", 



(35) 



where a s = g 2 /4irhc is the QCD strong coupling 
strength. Employing purely thermodynamic reasoning 
with the Gibbs-Duhcm equation, 



dP = sdT + rifd/if, 
f 

and the Euler relation for the energy density 



Ts-P + J2 



HfTlf, 



(36) 



(37) 



one finds the trace of the energy-pressure tensor 

9 eee u - 3P, 
B = Ts-4:P + ^2fi f n f , 



e = t 



dP 
df 



(38) 



If one replaces the quark chemical potentials with the 
quark activities 



z f = exp(/i//fc B T), 



(39) 
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FIG. 1: Shown are the relative pressure P/Pao and relative 
energy u/u x plotted as a function of temperature. 



then the thermodynamic identity for the trace of the 
energy-pressure tensor has the more simple form 



h rig) -ir. 



(40) 



or equivalently, the energy density is related to the pres- 
sure according to 



T 



dP 
df 



- P. 



(41) 



Consider the following simple model of a proposed 
quark-gluon plasma deconfinement. At very high tem- 
peratures, one considers the quarks and gluons to act 
as a gas exhibiting conformal symmetry. This implies a 
pressure P^ and an energy density obeying 



3P a 



const. 



(k B T) 4 



(hcf 



(42) 



At low temperatures this ideal behavior cannot persist 
and the conformal symmetry is broken by an energy gap 
A. The pressure may then be suppressed by a Boltzmann 
factor, 



P = Pcoe- A / fcBT , 



which, in virtue of Eq. (|41[) . implies 



-A/fc B T 



1 



A 



3k B T 



(43) 



(44) 



Plots of the quark-gluon deconfinement based on this 
model are shown in FIG. [TJ This simple model of the 
crossover between the confined phase and the quark- 
gluon plasma phase is in qualitative agreement with lat- 
tice QCD computations [23. The theoretical equilibrium 
statistical thermodynamic arguments presented here are 
not directly applicable to ultra high energy heavy ion 
nuclear collision for two reasons: (i) Very energetic colli- 
sions involve non-equilibrium processes and (ii) the col- 
lision fragments are most often described by QCD string 
fragmentation models [26j. 



VI. LUND STRING MODEL 

To compute non-equilibrium processes wherein quarks 
are connected to one another via strings (gluon electric 
flux tubes) we begin with the activation entropy Si(E) 
of a single string; It is@ 



Si_ 

k B 



E 



k B T 
S_ _ 1 
k B ~ 2 



E 



k B T Q 



S_ 

k B 



ln(3) 



(45) 



wherein the Hagedorn temperature To is related to the 
string tension a via 



k B T = 



3hca 

47T 



207 McV. 



(46) 



The numerical value of a (and thereby T ) is found from 
the experimental slope of the lowest meson Regge tra- 
jectory. The energy of a single string as a function of 
temperature may be computed via the microcanonical 
statistical thermodynamic law, 



1 
T 



dSjjE) 
dE : 



yielding, 



T 



T-Tn 



where $ 



7k B T 



(47) 



725 MeV. (48) 



The plot of energy as a function of temperature is given 
in FIG. [H Note that the energy always decreases as the 
temperature increases. The heat capacity is negative in 
the physically allowed string temperature range; i.e. 



dE 
C ~dT 



To 



< 0, 



LT-T _ 

> T >= —00 and 00 > T > T . 



(49) 
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FIG. 2: Shown is the energy of a single Lund string as a func- 
tion of temperature. Note that the heat capacity of the string 
C = (dE/dT) < in the physical but metastable temperature 
ranges -00 < T < and T < T < 00. 
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The physical interpretation of Eqs. (|49[) is worthy of note. 
Since in accordance with the second law of thermody- 
namics, stable systems have a positive heat capacity we 
thereby deduce that a single string during a collision is 
in a metastable state. 

As discussed in previous work[l 03, the string near 
the beginning of the collision has a very low energy 
< E < $. Thus, the initial string temperature is 
negative. The string cannot have negative energy which 
implies that string temperatures between zero and the 
Hagedorn temperature are not allowed; i.e. < T < Tq 
is unphysical. As the energy of the string energy increases 
due to the absorption of collision energy E — > $ — + , 
the temperature decreases T — > — oo. When only a lit- 
tle bit more collision energy is absorbed by the string 
E = $ + + the string temperature is virtually a positive 
infinite value. Any further cooling of the temperature 
requires more and more addition of the collision energy 
to be absorbed by a single string. We note in passing, 
that the notion of a negative temperature going positive 
by first approaching — oo and reappearing at +oo is not 
new. This property experimentally appears in laboratory 
nuclear magnetic resonant spin systems and laser media 
with inverted energy populations. 

For energy slightly above <3>, the string temperature is 
virtually infinite; i.e. T — > +oo as E — > $ + 0+. The 
question then arises (in the ultra high collision energy 
per nucleon limit) as to how much collision energy can 
be transferred to a single string in order to bring the 
temperature down from its very high temperature state. 
As will be made clear in the next Sec lVIIl the cooling 
time exponentially grows larger with decreasing temper- 
ature. Furthermore, the higher the collision energy per 
nucleon, the shorter the duration of a collision. When 
the duration of the collision is comparable to the cooling 
time for a given hot string temperature, the temperature 
stops getting smaller. The string temperatures stay high. 
While it is certainly possible to create many new strings 
with the ultra high collision energy, e.g. mesons, it will 
not be easy to cool the strings down from having a very 
high temperature. Let us now be more explicit regarding 
the time-scales in a stringy liquid droplet; i.e. the new 
high energy liquid droplet "pancake" model. 



VII. LOW VISCOSITY AND THE IDEAL FLUID 

The role of the entropy in determining transition rates 
follows from the the calculation rules of quantum me- 
chanical transition rates; In particular one averages over 
initial states and sums over final states. In general the 
number of quantum states O is related to the entropy S 
by the Boltzman formula 



S = fcfj In Q. 



(50) 



entropy difference 



Clt 



St 



(51) 



Transition rates may be phase space dominated by the 
exponential entropy factors. The dominant entropy 
arises from bosonic gluon excitations of the QCD string. 
The Fermi or Bose nature of physical particles depend 
only on the number, respectively odd or even, of quarks 
tied together (in a polymer-like fashion) by the gluon 
color electric flux string. The thermal relaxation times 
for very hot T — > +oo strings to cool down to finite tem- 
peratures strings depends on the exponential of an en- 
tropy difference 



T 

Toe 



exp 



5*i — Sx 



kB 



From Eqs.([i5|) and (@7|) one finds 



Si 
k B 

wherein 

"5*100 



To 



2 V T - Tn 



In 



T 



(T - T Q ) 



S\ c 



1 + bx 



2n 



In 



In 3. 



(52) 



(53) 



(54) 



Thus the cooling time scale from Eqs. ([52f and |53|) may 
be written as 



T = T a 



T-Tp 
T 



7/2 



exp 



k B (T-T ) 



(55) 



The viscosity implicit in QCD inspired string models 
is of a form which follows from Eqs.® and (|55|) ; It is 



n 



T-T 
T 



7/2 



exp 



4> 



AT -To) 



(56) 



To compute the attempt frequency r"^ 1 for a string to 
be created, for example between a quark anti quark pair, 
we need to consider the rate with which a quark in the 
"negative vacuum Dirac-Fermi sea" is excited into a pos- 
itive energy state leaving behind a quark hole. The phase 
space element along the axis of the string connecting the 
quark and the hole is 



d z N 



~ dpdx dp / dxdt\ 
~ 2irh ~ ~~dt V 277/7, J 



(57) 



wherein the rate of change of momentum of the quark is 
the force, i.e. the string tension a. Thus, the attempt 
frequency per unit time per unit length to excite a quark 
anti quark string, i.e. meson, is 



The ratio of the number of final states to the number of 
initial states is then determined by the exponential of an 



d 2 N 
dtdx 



a 
2nh' 



(58) 
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For string of length L, the energy is E 
attempt frequency 



aL yielding the 



1 



dN 
~dt 



E 

2nh 



MIT 



T 



T 



(59) 



wherein Eq. (|48|) has been invoked. 

From Eqs. (f56]) and follows a central result of this 
work; i.e. the high energy nuclear fluid viscosity 



V = Videal 



T-T a 
T 



9/2 



exp 



k B (T~T ) 



(60) 



with the ideal quantum viscosity determined as 
V « Videal if T > T , 



2irh 



< 



pc 
3> 



n, 

.2 



(61) 



Eqs. (f60|) and ([61]) for the viscosity of a stringy QCD fluid 
phase is quite similar to the viscosity of a polymer chain 
glass. As the temperature is lowered from above to To, 
the viscosity becomes extraordinarily high and the fluid 
behaves more as a glass then as an ideal fluid. For the 
ultra high energy collisions of heavy nuclei, the duration 
of the collision is sufficiently short that T ^ T). For such 
a range of hot string temperatures, the fluid viscosity 
is extraordinarily low, i.e. the viscosity is the quantum 
limited value r\ w rjideal- The QCD fluid phase is thereby 
ideal. 



VIII. CONCLUSION 



We have derived Kubo formula sum rules which al- 
lowed for the determination of the transverse viscosity in 
terms of the relaxation time of the pressure tensor. We 
then considered the quantized unit of circulation which 
contributes to turbulent eddy currents in the fluid. If the 
the quantum turbulent Feynman number T is to exceed 
classical turbulent Reynolds number 1Z, then there exists 
a quantum lower bound to viscosity which has been pre- 
viously conjecturedflll- Nuclear matter, as it appears in 
ultra high energy heavy nuclei scattering, was thought 
by some workers to have the formal conformal symme- 
try associated with deconfined glue at low baryon den- 
sity. Here, the string fragmentation model was employed 
wherein the quarks and anti quarks are connected to one 
another by strings whose electric flux tubes consist of 
glue. The notion of a low baryon density nuclear fluid 
is thereby defined as a very stringy gluon liquid with 
not so many quarks and anti quarks. This stringy liq- 
uid is not quite the same as a deconfined quark gluon 
plasma. The viscous properties of a QCD stringy liquid 
is as described by Eqs. (|60"f and (f6~lj) . The stringy liquid is 
closely analogous to a polymer glass. While the viscosity 
grows quite high as the glass (Hagedorn) temperature is 
approached from above, for temperatures well above the 
glass temperature (T ^> T ), the viscosity is quite small 
and quantum limited. These small viscosity values imply 
and almost ideal ultra high energy nuclear fluid. 
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